Estimation of the Local Density of States on a Quantum Computer 
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We report an efficient quantum algorithm for estimating the local density of states (LDOS) on a quantum 
computer. The LDOS describes the redistribution of energy levels of a quantum system under the influence of a 
perturbation. Sometimes known as the "strength function" from nuclear spectroscopy experiments, the shape of 
the LDOS is directly related to the survivial probability of unperturbed eigenstates, and has recently been related 
to the fidelity decay (or "Loschmidt echo") under imperfect motion-reversal. For quantum systems that can be 
simulated efficiently on a quantum computer, the LDOS estimation algorithm enables an exponential speed-up 
over direct classical computation. 
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A major motivation for the physical realization of quantum 
information processing is the idea, intimated by Feynman, that 
the dynamics of a wide class of complex quantum systems 
may be simulated efficiently by these techniques 1 1]. For a 
quantum system with Hilbert space size N, an efficient sim- 
ulation is one that requires only Polylog(iV) gates. This sit- 
uation should be contrasted with direct simulation on a clas- 
sical processor, which requires resources growing at least as 
N"^. However, complete measurement of the final state on a 
quantum processor requires 0{N'^ ) repetitions of the quantum 
simulation. Similarly, estimation of the eigenvalue spectrum 
of a quantum system admitting a Polylog(iV) circuit decompo- 
sition requires a phase-estimation circuit that grows as 0{N) 
ji"]. As a result there still remains the important problem of 
devising methods for the efficient readout of those characteris- 
tic properties that are of practical interest in the study of com- 
plex quantum systems. In this Letter we introduce an efficient 
quantum algorithm for estimating, to 1 /Polylog( A^) accuracy, 
the local density of states (LDOS), a quantity of central inter- 
est in the description of both many-body and complex few- 
body systems. We also determine the class of physical prob- 
lems for which the LDOS estimation algorithm provides an 
exponential speed-up over known classical algorithms given 
this finite accuracy. 

The LDOS describes the profile of an eigenstate of an un- 
perturbed quantum system over the eigenbasis of perturbed 
version of the same quantum system. In the context of many- 
body systems the LDOS was introduced to describe the ef- 
fect of strong two-particle interactions on the single particle 
(or single hole) eigenstates 1^015100]. More recently, the 
LDOS has been studied to characterize the effect of imperfec- 
tions (due to residual interactions between the qubits) in the 
operation of quantum computers This profile plays a 

fundamental role also in the analysis of system stability for 
few-body systems subject to a sudden perturbation 1 10|, such 
as the onset of an external field, and has been studied exten- 
sively in the context of quantum chaos and dynamical local- 
ization Quite generally the LDOS is related to the 
survival probability of the unperturbed eigenstate 



and there has been considerable recent effort to understand the 
conditions under which the LDOS width determines the rate 
of fidelity deca y un der imperfect motion-reversal ("Loschmidt 
echo") oEllilEl. 

A number of theoretical methods have been devised to char- 
acterize the LDOS for complex systems. These methods in- 
clude banded random matrix models mod- 
els of a single-level with constant couplings to a "picket- 
fence" spectrum |4, 20], and perturbative techniques with par- 
tial summations over diagrams to infinite order lljlll . Under 
inequivalent assumptions these approaches affirm a generic 
Breit-Wigner shape for the LDOS profile. 
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However, the extent to which these methods correctly describe 
any real system is generally not clear jT^, E^], and there- 
fore direct numerical analysis is usually necessary. It is woith 
stressing here that direct numerical computation of the LDOS 
requires the diagonalization of matrices of dimension N, and 
therefore demands resources that grow at least as N"^. Of 
course only coarse-grained information about the LDOS is of 
practical interest since one cannot even store the complete in- 
formation efficiently for large enough systems. However, for 
generic systems there is no known numerical procedure that 
can circumvent the need to manipulate the N x N matrix 
in order to extract even coarse information about its LDOS. 
In this Letter we report a quantum algorithm which enables 
estimation of the LDOS to l/Polylog(A^) accuracy with only 
Polylog(A^) resources. 

To specify the algorithm we represent the unperturbed 
quantum system by a unitary operator U, which may corre- 
spond either to a Floquet map, or to evolution under a time- 
independent Hamiltonian, 



U = exp{—iHoT) 



(2) 



We represent the perturbed quantum system by the unitary op- 
erator U{a), which we express in the form. 



U{a) = ex'p{-iSV)U, 



(3) 
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where S is some dimensionless parameter and V is a Hermi- 
tian perturbation operator. The variable a denotes an effective 
"perturbation strength" taking into account both the parameter 
5 and the size of the matrix elements of the perturbation, 



(4) 



where the average is taken only over directly coupled eigen- 
states. Let U\<j)j) = eyip{—i(j)j)\<pj), and U{a)\(j)k{o')) = 
exp (— i^fc(cr))|(/)fe(cr)) denote the eigenphases and eigenstates 
of the unperturbed and perturbed systems respectively. The 
LDOS for the j'th eigenstate of U is then, 



(./)) =^P(,^fc(a)|0,) <5(^-(</,fc(a) -</>,)), (5) 



where the transition probabilities. 



(6) 



are the basic quantities of interest. 
The coarse-grained distribution, 



(7) 



is a just the sum over the probabiUties for those perturbed 
eigenphases (pk{o') lying within a band A;. This band is cen- 
tered about angle 2t:1/M, with width A = 2n/M, and the 
integer I ranges from to M — 1. Similarly, an averaging over 
neighboring unperturbed eigenstates is often carried out to re- 
move the effects of atypical states. The combination of both 
operations yields the probability distribution. 



(8) 



where the normalization constant A^„i is just the number of 
unperturbed eigenphases in the angular range A^. In prac- 
tice one must choose M to be 0(log(A'^)) since otherwise the 
measured LDOS rj would contain an exponential amount of 
information and therefore could not be processed efficiently. 
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FIG. 1 : Circuit diagram for measuring the local density of states, 
consisting of two successive phase-estimation circuits on different 
operators. The diagonal line denotes a bundle of qubits and the thick 
vertical bar denotes a projective measurement of the quantum state in 
the computational basis. The upper register contains niq = log(M) 
qubits and the operations on the lower register are applied condition- 
ally m times, where the integer me [0, M — 1] is determined from 
the binary representation of the computational basis states in the up- 
per register. 



We now describe the algorithm for estimating the LDOS on 
a quantum processor. The circuit for this algorithm is depicted 
in Fig. 1. The lower register implements the perturbed and 
unperturbed maps U and U{(t), requiring Uq = 0(log2(iV)) 
qubits. The upper register holds the ruq = logj (M) ancillary 
qubits which fix the precision of the phase-estimation algo- 
rithm. The upper register always starts out in the 'ready' state 
|0). The appropriate choice of initial state p in the Uq register 
will depend on the context, as explained below. For the mo- 
ment we assume the lower register is prepared in a pure state, 
p = \'ijjo){i^o\- The first step of the algorithm involves esti- 
mating the eigenphases of the unperturbed operator U. This 
takes the initial state through the sequence. 



M-l 



|0) ® \^o) 



]= ^ |m)|^c,) 

m=0 

M-l 

= ^ \m){Ur\^o) 



m=0 

M-l Af-1 



= ~m E 1"^) E Cj (i:x.p{i(t>jm)\(t>j) 



m=0 j=0 



Af-1 
j=0 



(9) 



where c, 



, IV'o)- The state rrij is the nearest m^-bit binary 



approximation to the j'th eigenphase of U 



(10) 



Upon strong measurement of the register one obtains and 
records a single outcome m, and the state of rig register must 
then be described by (viz, 'collapsed to') the updated pure 
state. 



(11) 



corresponding to the subspace of eigenstates with eigenphases 
in the band A^ of width A = 27r/M about the phase 
2-Km/M. To keep normalization the coefficients have been 
rescaled as follows. 



(12) 
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Next we reset the qubit register to the ready state and run 
the phase-estimation algorithm on the operator U {a), produc- 
ing the final state. 



AT-l 



IV') = E E ^(^'•^') i*^*^) \<i>k{(^)), 



(13) 



where (f>k{l) = 2TTmk/M is an m^-bit approximation to 
(/)fc(cr). The complex coefficients b{k\j) = {4>j\(t>k{'^)) are 
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the inner product of the perturbed and unperturbed eigen- 
states. Measurement of the niq register now reveals an out- 
come /, associated with the eigenphases in the angular range 
2ttI/M ± a/2. The outcome / occurs with probability, 



(14) 



which is conditional on the earlier outcome m and the choice 
of initial state. 

We now specify how the initial state may be chosen to elim- 
inate unwanted fluctuations arising from the variables Cj in 
Eq. . Before describing the general solution we consider 
first a special case of particular interest: when a known eigen- 
state of U may be prepared efficiently. Such an initial state 
may be prepared (or well approximated) by an efficient cir- 
cuit when U consists of some sufficiently simple integrable 
system (e.g., a non-interacting many-body system). In this 
case we have cj — Sjk, and the final probability distribution 
Eq.^^reduces exactly to the (coarse-grained) kernel Eq.0 



P^,(Z|m)->P(A, 



(15) 



When the eigenphase associated to the prepared eigenstate is 
known to sufficient accuracy (so that m is known), it is not 
even necessary to perform the first phase estimation routine. 
In the general case of a generic quantum system, it is sufficient 
to prepare the maximally mixed state as the initial state, in 
which case the final probability distribution reduces exactly 
to the (coarse-grained and averaged) probability kernel Eq.|8l 
i.e.. 



Pi/jv(Z|m) 



1 
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(16) 



This probability kernel contains all the information needed 
to compute the (coarse-grained and averaged) LDOS, 
Tjmi^Trk/M) = P(/|m)(5fc_(;_,„), completing our deriva- 
tion. 

The algorithm described above remains efficient provided 
that the quantum maps U and Us admit Polylog(A^) gate de- 
compositions. Such decompositions have been identified both 
for many-body systems with local interactions and for a wide 
class of few-body quantized classical models. As mentioned 
earlier, for practical purpose M should be Polylog(Af) so the 
overall circuit of Fig.^is indeed efficient for such systems. 

We now turn to the question of how many times K the al- 
gorithm must be repeated to arrive at interesting physical con- 
clusions about the final probability distribution. This issue 
arise because the final probability distribution is not measured 
directly on the quantum processor; rather, it governs the rela- 
tive frequency of outcomes obtained in each repetition of the 
algorithm. Indeed, it is by repeating the algorithm illustrated 
in Fig.^and accumulating joint statistics of the I and m out- 
puts that one can estimate the parent distribution P(/|m). The 
accuracy of this estimation depends on the number of times 



K the distribution is sampled. In order to bound K it is con- 
venient to cast the physical problems related to the LDOS in 
terms of hypothesis testing. We consider the important case of 
testing which of two candidates distributions 771 or 772 best de- 
scribes the LDOS of a given system and a given perturbation. 
For example, one might be testing whether the Lorentzian 
has one of two candidate widths, or whether the profile is 
Gaussian or Lorentzian. Only when K < Polylog(A^) will 
the overall computation remain efficient. This problem is re- 
solved in general by the Chernoff bound [22]. A random vari- 
able is distributed according to either Pi{x) or P2{x), and we 
wish to determine which distribution is the right one. Then, 
the probability P^ that we make an incorrect inference de- 
creases exponentially with the number of times K the variable 
was sampled: Pg < . Here, < A < 1 is a measure of 
similarity between distributions defined as 

A== min VPi(x)"P2(a;)(i-"); (17) 

X 

in particular, A is bounded above by the fidelity between 
Pi and P2. Thus, a constant error probability e requires a 
sample of size K — log(e)/ log(A). Therefore, as long as 
the concerned distributions are at a Polylog(iV) distance, i.e. 
1 — A > 1 /Polylog(A^), they can be distinguished efficiently. 
We note that the test can be inconclusive when both hypothe- 
sis are equally likely to describe the underlying physics. 

Efficient application of the LDOS algorithm under these 
restrictions may be illustrated explicitly by working through 
a problem of practical interest from the recent literature. We 
consider the problem of testing whether the Breit-Wigner pro- 
file Eq.Q applies when a given quantized classically chaotic 
model is subjected to a perturbation of interest. From the BGS 
conjecture \2m and studies of (banded) random matrix models 
II3L I171II8LI19I1 . it is generally expected that for fully chaotic 
models with generic perturbations Eq.^applies with. 



r(cr) = 2TTa^pE, 



(18) 



provided that the effective perturbation strength hes in the 



range. 



1 < crpE < Vb, 



(19) 



where b is the bandwidth of the perturbation in the ordered 
eigenbasis of U and is the level density. It should be 
stressed that F may be estimated a priori if the perturbation is 
known I.14».. 17^ .181. Deviations from this hypothesis can arise 
for a wide variety of reasons (i.e., integrable or mixed clas- 
sical dynamics in the unperturbed or perturbed system, non- 
generic properties of the perturbation, hidden symmetries, etc) 
and therefore analysis of the LDOS remains an active area of 
numerical study for both dynamical models and real 

systems |5]. 

The lower bound of Eq. ^]is determined from the break- 
down of perturbation theory and leads to a width F that de- 
creases linearly with A^. Since the circuit can only efficiently 
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resolve the LDOS with accuracy l/Polylog(A'^), the BW pro- 
file with width F ~ 0{N^^) may not be verified efficiently 
near this lower bound. However, near the upper bound of 
Eq.[^the validity of the BW profile may be tested efficiently. 
In the case of fully chaotic models one has b — N/2 and 
the upper bound for F is therefore 0(1). Hence the validity of 
Eq.fflprovides a hypothesis which may be tested efficiently for 
any perturbation such that l/Polylog(iV) < F(cr) < 0(1). 
Near this bound one can also determine whether the chaotic 
model exhibits dynamical localization, since in this case one 
has a bandwidth b ^ N/2 and the LDOS will cease to main- 
tain the BW profile when bp]^^ < F((t) < 0(1). Indeed for 
some models the localization length / of the eigenstates scales 
as / ~ lEsll . and hence this length may be estimated 

using the LDOS algorithm with only Polylog(A^) resources. 

In summary we have reported an algorithm for efficiently 
estimating the LDOS of a quantum system subject to pertur- 
bation. There is wide range of contexts in which important 
coarse features of the LDOS, such as the width, may be es- 
timated with only Polylog(A^) resources. We have described 
in detail the important problem of testing the Breit-Wigner 
hypothesis as one example for which the LDOS estimation al- 
gorithm gives an effective exponential speed up over classical 
computation. 
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